Symmetry classification of energy bands in graphene and silicene 
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We present the results ol the symmetry classification of the electron energy bands in graphene and 
silicene using group theory algebra and the tight-binding approximation. The analysis is performed 
both in the absence and in the presence of the spin-orbit coupling. We also discuss the bands 
merging in the Brillouin zone symmetry points and the conditions for the latter to become Dirac 
points. 

PACS numbers: 73.22.Pr 



INTRODUCTION 

Since graphene was first isolated experimentally [l|, 
it is in the focus of attention of both theorists and ex- 
perimentalists. Obviously, understanding of the symme- 
tries of the electrons dispersion law in graphene is of cru- 
cial importance. Actually, the symmetry classification 
of the energy bands in graphene (or "two-dimensional 
graphite") was presented nearly 60 years ago by Lomer 
in his seminal paper [3|. Later the subject was analyzed 
by Slonczewski and Weiss [3], Dresselhaus and Dressel- 
haus [J], Bassani and Parravicini [5|. Recent approaches 
to the problem are presented in the papers by Malard et 
al. [6| , Manes [7| and in our publication [8| . 

The present work has two aspects: a pragmatic and a 
pedagogical one. The first aspect is connected with the 
recent synthesis of silicene, the counterpart of graphene 
for silicon, with buckled honeycomb geometry. This novel 
two-dimensional material has attracted considerable at- 
tention both theoretically [9|, llO| and experimentally re- 
cently, due to its exotic electronic structure and promis- 
ing applications in nanoelectronics as well as its compati- 
bility with current silicon-based electronic technology. So 
we present the symmetrty analysis of the silicene electron 
bands. 

The pedagogical aspect is connected with the fact that 
different approaches to the symmetry classification, even 
if giving the same results, are based on different meth- 
ods of applications of group theory. Thus in our previous 
paper [8] the labeling of the bands was based on com- 
patibility relations and guesses. In the present work we 
show that in the framework of the tight-binding approxi- 
mation the representations of the little group in the sym- 
metry points can be rigorously found in the framework 
of the group theory algebra. Though the idea of using 
the tight-binding approximation is by no means new (it 
was used already in the work by Lomer), our mathemat- 
ical approach is totally different, as one can easily see 
comparing the present work with [4], and, to our opin- 
ion, more convenient for applications. This statement is 
supported by the analysis of the symmetry of the energy 
band in silicene. 

We also generalize the symmetry classification by 



taking into account the spin-orbit coupling both for 
graphene and for silicene. This, to the best of our knowl- 
edge, wasn't done before even for graphene. Though in 
graphene the spin-orbit coupling is very weak, the prob- 
lem is interesting in principle. One can expect that in 
silicenr the coupling is stronger, and it will become even 
more so for graphene related materials from heavier ele- 
ments, provided they can be synthesized. 

To remind to a reader a few basic things, important for 
the symmetry classification of the bands in any crystal, 
consider a point sub-group R of the space group charac- 
terizing the symmetry of a crystal (we restrict ourselves 
with the consideration of symmorphic space groups). 
Any operation of the group R (save the unit transfor- 
mation) takes a general wavevector k into a distinct one. 
However, for some special choices of k some of the oper- 
ations of the group R will take k into itself rather than 
into a distinct wavevector. These particular operations 
are called the group of k; it is a subgroup of the group R. 
Points (lines) in the Brillouin zone for which the group of 
the wavevector contains elements other than the unit el- 
ement are called symmetry points (lines). We may use a 
state (states) corresponding to such a special wavevector 
to generate a representation for the group of k [llj, Il2| . 
In this paper we consider crystals with the hexagonal 
Brillouin zone. In this case the symmetry points are F - 
the center of the Brillouin zone, the points K which are 
corners of the Brillouin zone and the points M which are 
the centers of the edges of the Brillouin zone. 



TIGHT-BINDING MODEL 

We'll deal with the materials with a basis of two 
atoms per unit cell, and we'll search for the solution 
of Schrocdinger equation as a linear combination of the 
functions 



V-^' 



/3;k 



Ev^^(^ 



R.T 



^ik-R,, 



(1) 






where '0/3 are atomic orbitals, j — A, B labels the sub- 
lattices, and Rj is the radius vector of an atom in the sub- 
lattice j. A point symmetry transformation of the func- 



tions tpg.-i^ is a direct product of two transformations: the 
transformation of the sub- lattice functions 



ij^' , where 
(2) 



R, 



and the transformation of the orbitals ipp. Thus the rep- 
resentations reahzed by the functions ([1]) will be the di- 
rect product of two representations. Generally, this rep- 
resentation will be reducible. To decompose a reducible 
representation into the irreducible ones, it is convenient 
to use equation 



^E^(G)x:(G), 



(3) 



G 



which shows how many times a given irreducible repre- 
sentation is contained in a reducible one [11]. Additional 
information about the representations can be obtained if 
we use projection operator [14] 



Oo = yEx:(G)P(G), 



G 



(4) 



where Ha is the dimensionality of the irreducible repre- 
sentation a and P is the operator corresponding to a 
given transformation G The operator projects a given 
function to the linear space of the representation a. For 
a one dimensional representation the operator thus gives 
basis of the representation. 



GROUP THEORY ANALYSIS IN THE 

TIGHT-BINDING MODEL WITHOUT THE 

SPIN-ORBIT COUPLING 

Our tight-binding model space includes four atomic or- 
bitals: |s,p >. (Notice that we assume only symmetry 
of the basis functions with respect to rotations and re- 
flections; the question how these functions are connected 
with the atomic functions of the isolated carbon atom is 
irrelevant.) 



Graphene 

The Hamiltonian of graphene being symmetric with 
respect to reflection in the graphene plane, the bands 
built from the \z > orbitals decouple from those built 
from the \s,x,y > orbitals. The former are odd with 
respect to reflection, the latter are even. In other words, 
the former form tt bands, and the latter form a bands. 

The group of wave vector k at the point F is D(ih, at 
the point K is D^h , at the lines F — isT is C2v [Sl, ll5| . 
The group of wave vector k at the point F is Deh , a t 
the point K is D^h , at the lines T — K is C2v la |l5| . 
The representations of the groups D^h and Dqii can be 
obtained on the basis of identities 
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E 2C3 3U2 
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1 1 1 
1 1 -1 
2-10 



C(iv 


E C2 2C3 26*6 3cr„ So-; 


Ai 


11 1 111 


A2 


11 1 1-1-1 


B2 


1-1 1-1 1-1 


Bi 


1-1 1-1-1 1 


E2 


2 2-1-1 


El 


2-2-1 1 



Table I: Character table for irreducible representations 
of -D3 and Cqv point groups 



The irreducible representations of the groups Dy, and Cgi, 
are presented in the Table 1. Each representation of the 
group, say Ai, begets two representations of the group 
Z?3/i : A'l and A'l] prime means that the representation 
is even with respect to reflection a in double prime means 
that it is odd. 

The irreducible representations of the group Cg^ are 
also presented in the Table 1. Because the inversion 
transformation / can be presented as 



/ = C20'/i, 



(6) 



the representations of the group can be classified as sym- 
metric (g) or antisymmetric (u) with respect to inversion. 
Thus each representation of the group Cq^ , say Ai , begets 
two representations: Aig and Ai„. 

Notice that the orbitals |s > (or |z > ) realize Ai 
representation both of the group D^ and of the group 
Cgi; , hence the representations of the groups realized by 

A R 

the functions tp^ ^.j. will be identical to those realized by 

the sub-lattice functions (j)^^' . 

Let us start from the symmetry analysis at the point 
F. Because the transformations C2,Ce,av change sub- 
lattices, the characters corresponding to these transfor- 
mations are equal to zero. The transformations -B, C3, ay/ 
leave the sub-lattices as they were. Hence from Table U 



we see that the functions 
sentation 



iA,B 



realizes reducible repre- 



i?r - Ai + B2 



(7) 
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Cfii, X C's 



(5) 



of the group Cev ■ 

Taking into account the symmetry of the states relative 
to reflection in the plane of graphene , we obtain that 
at the point F the functions ipz (here and further on, 
when this is not supposed to lead to a misunderstanding, 
we'll suppres the index k in ipp-u) realize Aiu and i?2g 
representations of the group Dgh, characterizing tt bands, 
and the functions ips realize Aig and A2U representations 
of the group, characterizing a bands. 

Acting by projection operators Oa and Ob on a func- 
tion ip^ , we obtain that the irreducible representation Ai 
is realized by symmetric combination of the A and B 
orbitals, and the irreducible representation B - by the 
antisymmetric combination. One can expect that the 
flrst case occurs in the hole band, and the second - in the 
electron band. 



The orbitals \x,y > realize representation E-^ of the 
group Cqv [13| . Hence, representation of the group real- 
ized by the functions V'^'^ can be decomposed as 



Ei-kRt = Ei+ E2 



(8) 



Taking into account the symmetry of the states relative 
to reflection in the plane of graphene, we obtain that 
at the point T the functions '4'x,y realize Eig and £'211 
representations of the group D^h, characterizing a bands. 
To find wavefunctions realizing each of the irreducible 
representations we apply the projection operators and 
obtain 



Ofi.V'i 



(V'i 



>i), 



Oe,^'± 



(V'i 



V'i) (9) 



where j = B ii j = a, and vice versa. Thus represen- 
tation Ei{E2) is realized by symmetric (antisymmetric) 
combinations of \x > and \y > orbitals. One can expect 
that the first representation is realized at the hole band, 
and the second - at the valence band. 

Now let us perform the symmetry analysis at the point 
K. The orbitals \s > (or I2; > ) realize Ai representation 
of the group D^^. 

The representation of the group realized by the func- 
tions (/)jj.' is determined by the transformation law of 
the exponentials e*^'^^ under the symmetry operations. 
Rotation of the radius vector by the angle 27r/3 anti- 
clockwise, is equivalent to rotation of the vector \bfK 
in the opposite direction, that is to substitution of the 
three equivalent corners of the Brillouin zone: K — > 
K2 ^ K3 ^ K where K = (27r/3a, 27r/3A/3a) , K2 = 
(0, -47r/3V3a) and K3 = (-27r/3a, 27r/3V3a). The ro- 
tation multiplies each basis vector by the factor e^^*/'^. 
Using Eq. [3]) we obtain 



as = if2-e2W3 



-27ri/3^ 



1, 



(10) 



Hence the functions (p-^ realize irreducible representa- 
tion E of the group -D3. Taking into account the sym- 
metry of the states relative to reflection in the plane of 
graphene, we obtain representation E" of the group D^h, 
realized by the tpz functions, characterizing tt bands, and 
representation E' , realized by the V's functions, charac- 
terizing a bands. 

Consider now representation of the group D^ realized 
by the quartet of functions ■i/'^.y at the point K. The or- 
bitals \x,y > realize representation E of the group [13J . 
Hence, representation of the group realized by the func- 



tions tp^'if can be decomposed as 



Ex E = A1+A2+E. 



(11) 



Taking into account the symmetry of the states relative 
to reflection in the plane of graphene, we obtain repre- 
sentations A'^,A'2 and E' of the group D^h, realized by 
the ipx.y functions, characterizing a bands. 



Civ 


D(.h Dsh 


Rep 


Compatible with 


A2 
Bi 
B2 


Alg,B2u, Elu, E2g 
Alu,B2g, Elg,E2u 
Blu, A2g, Elu, E2g 

A2U , Big, Elg, E2u 


A'i,E' 
A'i,E" 
A'2,E' 
A'i,E" 



Table II: Compatibility relations 



Acting by projection operators, we obtain that the rep- 
resentation Ai is realized by the vector space with the 
basis vector ■0^ -I- V':? , and the representation A2 is real- 
ized by the vector space with the basis vector ^^ — tp^. 
The vector spaces realizing representations Ai and A2 be- 
ing found, the representation is obviously realized by the 
vector space spanned by the vectors V'+ — "0- ; V"- + '/'+ ■ 

Because the irreducible representation E' is realized 
both by ips and ipx,y functions, these representations 
should be considered together. According to Wigner 
theorem [16| we still have two representations, each of 
them being realized by two functions from a quartet 
i/'^'^ , "0^ — TpE , ip- + "0+ • Each E' representation charac- 
terizes two a bands, merging at the point K. The sym- 
metry of the electron bands at the points F and K being 
determined, the symmetry at the lines T ~ K follows un- 
equivocally from the compatibility relations, presented 
in Table [8J, ll5|. The table shows compatibility of the 
representations of the point group , realized at the sym- 
metry line , with those realized at the symmetry points 
and . The results of this Section are presented on Fig. (TJ 
reproduced from [8|. 
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Figure 1: (Color online) Graphene band structure 

evaluated with use of the FP-LAPW method. The 

dashed line shows the Fermi energy. [Reproduced from 

Ref. i.l 



Rep 


E 2C3 3(/2 / 2S6 ?>(Ja 


Ai„ 


111111 


A2c 


1 1-11 1-1 


Ai„ 


1 1 1-1-1-1 


A2u 


1 1-1-1-1 1 


E, 


2-1 2-1 


Eu 


2-1 0-2 1 



Table III: Characters tabic for irreducible 
representations of D^d point group 



Silicene 

The difference between silicene (or symmetrically 
equivalent to it buckled graphene) and graphene for our 
consideration is due solely to the decreased symmetry of 
the former. The group of the wavevector at the point 
K in silicene is D3, at the point F - D-^d (this is also the 
point group of silicene). The representations of the group 
D^d we can obtain on the basis of identity 



D:id = ^3 X C,. 



(12) 



The direct product has twice as many representations as 
the group Da, half of them being symmetric (denoted by 
the sufhx g), and the other half antisymmetric (sufhx u) 
with respect to inversion. The characters of the repre- 
sentations of the group are presented in the Table IIIII 

The symmetry analysis in silicene parallels that in 
graphene, so we'll be brief. 

At the point K the functions (p-^ realizes irreducible 
representation E of the group D^. 

Orbitals \s > and \z > realize Ai representation, and 
the orbitals \x,y > realize E representation of the group 
D3. Thus at the point K the functions "03 realize repre- 
sentation E of the group Z^a, and the same same can be 
said about the functions tpz- Reducible representation 
realized by the functions ipx.y can be decomposed into 
the irreducible ones: 



Ex E = A1+A2+E. 



(13) 



So when the symmetry is reduced by going from graphene 
to silicene, the representations A'^ and A2 turn into Ai 
and A2. Representation E" and two representations E' 
turn into three representations E. Loosing the reflection 
in plane symmetry, we can not claim now that one repre- 
sentation is realized exclusively by |z > orbitals. All the 
E representations mix |s,p > orbitals. 

At the point F the functions 
representation of the group D3: 



Q ' realizes reducible 



Rr = Aiq + A 



2m- 



(14) 



Orbitals \s > realize Aiu, orbitals \z > - A2U, and or- 
bitals \x,y > - Eu representations of the group. Thus 



the wavefunctions m realize reducible representation of 
the group D^d which is 



(Al„ + A2u + Eu) X {Aig + A2u) 

= Alu + A2g + A2u + Aig + Eg + Eu 



(15) 



So when the symmetry is reduced by going from graphene 
to silicene, the representations E2g and Eiu turn into the 
representations Eg and E^ respectively. 

The band structure of silicene is being different from 
that of graphene, the merging of the bands is no different. 
The statement becomes clear when comparing Fig. [21 
reproduced from Ref. [10| , with Fig. [TJ 




Figure 2: (Color onhne) DFT-PBE band structures for 

silicene. The dashed line shows the Fermi energy and 

the insets show the spectrum near the Fermi level in the 

vicinity of the K point. [Reproduced from Ref. [lO| 

with permission.] 



GROUP THEORY ANALYSIS IN THE 

TIGHT-BINDING MODEL WITH THE 

SPIN-ORBIT COUPLING 

In the absence of spin-orbit coupling, electron spin can 
be taken into account in a trivial way: each band we 
considered was doubly spin degenerate. When the spin- 
orbit coupling is taken into account, the algebra, actually, 
is not that different. The symmetry, and the represen- 
tations realized by the sub-lattice functions ^ remaim 
the same. However, instead of atomic orbitals we should 
consider atomic terms: /3 enumerates states terms from 
doublets \s >(^/^^ \p >(i/2) and quartet \p >(3/2), 

Due to the semi-integer value of J we have to consider 
double-valued representations realized by the atomic 
terms (and by the crystal wave functions). To find (the 
characters of) two- valued irreducible representations, it 
is convenient to introduce the concept of a new element 
of the group (denoted by Q); this is a rotation through an 



D',u 


E Q an C3 C'i S3 S'i 3U2 3(7„ 
anQ ClQ CiQ SlQ S3Q W2Q 3a„Q 


T7 


2-2 1-1 -v/3 V3 
2-2 1 -1 x/3 -V3 
2-20-22 



Table IV: Characters for two-valued irreducible 
representations of group Dg^ 



Rep 


E Q C3 C'i 3U2 SU2Q 
CiQ CzQ 


r4 


1-1-1 1 i -i 
1-1-1 1 i -i 
2-2 1-10 



Table V: Characters for two-valued irreducible 
representations of group D'^ 



angle 2tt about an arbitrary axis, and is not the unit el- 
ement, but gives the latter when applied twice: Q^ — E. 
The characters of the rotation by angle applied to 
the term I . . . >'^-') is 



xM) 



sin(J+^] 



With respect to the inversion / the character is 

X,(/)-±(2J+l), 



(16) 



(17) 



where the sign plus corresponds to the s states, and the 
sign minus to the p states. Finally, the characters corre- 
sponding to reflection in a plane a and rotary reflection 
through an angle (j) are found writing these symmetry 
transformations as 



Silicene 

The two- valued representations of are presented in Ta- 
ble|V][il. 

Each of the doublets \s >^^/'^\\p >(i/2) realizes r4 
representation of the group. Quartet \p >(^/2) realizes 
this representation twice. The sub-lattice functions real- 
ize representation E of the group. From Equation ([3]) we 
obtain 

r4xs = r4 + r5 + r6. (20) 

(For the same reasons as for ordinary representations, 
two complex conjugate two-valued representations must 
be regarded as one physically irreducible representation 
of twice the dimension.) Thus at the point K four bands 
realize representation r4 of the group D'^ each, and four 
bands realize representation r5 4- Fg each. 



= /C2 



5(0) = /C((/. + ^). 



(18) 



Both in graphene and in silicene we'll restrict ourselves 
by the symmetry analysis at the point K. 



Graphene 



The sub-lattice functions 



,A.B 



realize Eg representa- 



tion of point group D^h- The electron terms realize two- 
valued representations of the group, which are presented 
in Table [IV] [13. 

Doublet \s >(i/2) realizes T-j representation of the 
group; doublet \p >(^"-' realizes Fg representation, quar- 
tet ^p(3/2) realizes Fy representation twice. (We decided 
to use chemical notation for the single-valued represen- 
tation, and BSW notation for double-valued representa- 
tions [ 11]].) The sub-lattice functions realize representa- 
tion Eg of the group; from Eq. ^ we obtain 



Fr X £;„ = Fg X ^„ = F7 + Fg 



(19) 



Thus at the point K four bands realize representation F7 
of the group each, and four bands realize representation 
Fg each. In particular, we obtained the (well known) 
result that the four-fold degeneracy (including spin) of 
the bands merging at the point K is partially removed 
by the spin-orbit coupling, and only two-fold (Kramers) 
degeneracy is left. 



DIRAC POINTS 

In this final part of the paper we would like to clarify 
the relation between the symmetry and the existence of 
Dirac points. According to the classical approach [ISlJigj . 
the merging of the bands at a point ko is connected with 
the multi (higher than one) - dimensional representation 
of the space group Go, realized in this point. Looking for 
a linear dispersion point in the vicinity of the merging 
point we may use the degenerate k • p perturbation the- 
ory. Let a two-dimensional irreducible representation is 
realized at a point kg. Expanding the wavefunction with 
respect to the basis of the two-dimensional irreducible 
representation 



V'(k)=^c,(k)V^,(ko), 
for the expansion coefficients we obtain equation 

^^ilP^c,(k)=e(k)c,(k), 

where p^ =< V'i(ko)lpl^j(ko) > (of course, we need the 
absence of inversion symmetry at the point, for the ma- 
trix elements to be different from zero). The dispersion 
law is given by the equation 



(21) 



(22) 



e(k) 



/ ^ Ojci'^a i 



I'^lajskakp, 



(23) 



ap 



where a,/3 are cartesian indexes x,y. Eq. ([23)1 should 
contain only combinations of wavevector components 
which are invariant with respect to all elements of the 
group Gq. In the case when the group Gq does not have 
any vector invariants, and the only tensor invariant is the 
quantity k^ + k^, we obtain the dispersion law 



e(k) = ±vk, 



(24) 



which, like it was shown by Dirac himself in 1928, guar- 
anties that Eq. (P^ is Dirac equation, in the sense the 
the matrices p^. and Py satisfy anticommutation relations 



pI^pI'^ I, 



{Px,Py}+ = 0, 



where / is the unity matrix. 

To be more specific, consider the groups of wavevector 
at the point F; in graphene it is Ds^, and in silicene it is 
Z?3. In both cases, to find the dispersion law at the point 
r it is enough to study invariants of the group D^. And 
we can easily check up that both conditions, necessary 
for the existence of the Dirac point, are satisfied. And 
we can easily check up that both conditions, necessary 
for the existence of the Dirac point, are satisfied. This 
explains, in particular, why the band calculations show 
the existence of Dirac points in silicene [9|, llO| , which has 
lower symmetry than graphene. 

In general, the role of the tight binding approximation 
in symmetry classification of the bands in graphene, like 
its role in symmetry classification of bands in other crys- 
tals, is only auxiliary. The approximation greatly helps in 
the classification and sheds additional light on the nature 
of the bands, but one must remember that that there are 
more important things that this or that approximation - 
and that is symmetry. 



CONCLUSIONS 

This paper presents an applications of group theory to 
very important cases of graphene and silicene. 



* [Eugene. Kogan@biu.ac.il| 
[1] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
Y. Zhang, S. V. Dubonos, I. V. Grigorieva, A. A. Firsov, 
Science 306, 5696 (2004). 
[2] W. M. Lomer, Proc. Roy. See. A 227, 330 (1955). 
[3] J. C. Slonczewski and P. R. Weiss, Phys. Rev. 109, 272 

(1958). 
[4] G. Dresselhaus and M. S. Dresselliaus, Phys. Rev. 140, 

A401 (1965). 
[5] F. Bassani and G. Pastori Parravicini, Nuovo Cim. B 50, 
95 (1967). 
(25) [6] L. M. Malard, M. H. D. Guimaraes, D. L. Mafra, M. S. C. 

Mazzoni, and A. Jorio, Phys. Rev. B 79, 125426 (2009). 
[7] J. L. Manes, Phys. Rev. B 85, 155118 (2012). 
[8] E. Kogan and V. U. Nazarov, Phys. Rev. B 85, 115418 

(2012). 
[9] C.-C. Liu, Hua Jiang, Yugui Yao, Phys. Rev. B84, 
195430 (2011). 
[10] N. D. Drummond, V. Zolyomi, and V. I. Falko, Phys. 

Rev. B85, 075423 (2012). 
[11] C. Kittel, Quantum Theory of Solids, (John Wiley & and 

Sons. Inc. New- York, London 1963). 
[12] W. A. Harrison, Solid State Theory, (McGraw Hill Book 

Company, New York-London-Toronto, 1970). 
[13] L. D. Landau and E. M. Lifshitz, Quantum Mechanics, 

(Pergamon Press, 1991). 
[14] R. S. Knox and S. Gold, Symmmetry in the Solid State, 

(W. A. Benjamin, Inc., 1964, New York, Amsterdam). 
[15] C. Thomsen, S. Reich, J. Maultzsch, Carbon Nanotubes: 
Basic Concepts and Physical Properties, (Wiley Online 
Library, 2004 WILEY- VCH Verlag GmbH). 
[16] M. I. Petrashen, E. D. Trifonov, J. L. Trifonov, Appli- 
cations of Group Theory in Quantum Mechanics (Dover 
Publications, 2009). 
[17] G.F. Koster, J.O. Dimmock, R.G. Wheeler, and H. Statz, 
Properties of the Thirty-Two Point Groups (MIT Press, 
Cambridge 1964). 
[18] F. Hund, Z. Physik 99, 119 (1936). 

[19] C. Herring, Phys. Rev. B52, 365 (1937); J. Franklin Inst. 
223, 525 (1942). 



ACKNOWLEDGMENTS 



Discussions with J. L. Manes, V. Falko and Hua Jiang 
were very illuminating for the author. 



